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ABSTRACT 

We study the stability of compressible cylindrical differentially rotating flow 
in the presence of the magnetic field, and show that compressibility alters quali- 
tatively the stability properties of flows. Apart from the well-known magnetoro- 
tational instability that can occur even in incompressible flow, there exist a new 
instability caused by compressibility. The necessary condition of the newly found 
instability can easily be satisfied in various flows in laboratory and astrophysical 
conditions and reads BgB^Q' ^ where B s and B^ are the radial and azimuthal 
magnetic fields, Q' = dfl/ds with s being the cylindrical radius. Contrary to the 
magnetorotational instability that occurs only if Q decreases with s, the newly 
found instability operates at any sign of Q'. The considered instability can arise 
even in a very strong magnetic field that suppresses the magnetorotational in- 
stability. 

Subject headings: magnetohydrodynamics - instabilities - accretion, accretion 
discs - galaxies: magnetic fields 



1. Introduction 

Magnetohydrodynamic instabilities and turbulence generated by these instabilities can 
play an important role in enhancing transport processes in various astrophysical bodies, 
such as accretion and protoplanetary disks, galaxies, stellar radiative zones, etc. Anomalous 
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turbulent transport caused by instabilities can be particularly important in differentially 
rotating and magnetized gas where a number of MHD instabilities may occur. It is known 
since the classical papers by Velikhov (1959) and Chandrasekhar (1960) that a differentially 
rotating flow with a negative angular velocity gradient and a weak magnetic field is unstable 
to the magnetorotational instability. This instability has been analyzed in detail for stellar 
conditions (see Fricke 1969, Acheson 1978, 1979). In the context of accretion disks, this 
instability was first considered by Balbus & Hawley (1991). There are strong arguments 
that the necessary condition of the magnetorotational instability (a decrease of the angular 
velocity with cylindrical radius) is fulfilled in disks (see, e.g., Kaisig, Tajima, & Lovelace 1992, 
Kumar, Coleman, & Kley 1994, Zhang, Diamond, & Vishniac 1994). The instability exists 
not only for short wavelength perturbations, but also for global modes with scales comparable 
to the disk height (Curry, Pudritz, & Sutherland 1994, Curry & Pudritz 1995). Likely, the 
magnetorotational instability may develop in astrophysical plasma under a wide spectrum of 
conditions. The effect of radiation on stability of magnetized differentially rotating flow has 
been studied analytically by Blaes & Socrates (2001) and numerically by Turner, Stone, & 
Sano (2002) who argued that the stability criterion remains unchanged in this case but the 
growth rate can be lower. Recently, stability properties of magnetized accretion flows have 
been considered in the Kerr metric by De Villiers & Hawley (2003), De Villiers, Hawley, & 
Krolik (2003), and Hirose et al. (2004) who found that the magnetorotational instability can 
drive turbulence in this case as well. 

The numerical study of the magnetorotational instability in magnetized accretion disks 
both in local (Brandenburg et al. 1995, Hawley, Gammie & Balbus 1995, Matsumoto & 
Tajima 1995, Arlt & Riidiger 2001) and global (Armitage 1998, Hawley 2000) approximations 
show a breakdown of laminar Keplerian flow into well-developed turbulence. The generated 
turbulence can be an efficient mechanism of the a- viscosity and angular momentum transport 
in accretion disks (Brandenburg et al. 1996, Hawley, Gammie, & Balbus 1996, Stone et al. 
1996). At low temperatures and high densities, however, the level of turbulence induced 
by the magnetorotational instability can decrease significantly or even turns off completely 
(Gammie 1996, Gammie & Menou 1998). 

In recent years, many simulations of accretion disks have been performed, and much of 
the dynamics was interpreted as a direct consequence of the magnetorotational instability. 
Obviously, however, that this instability cannot be the only instability that operates in 
differentially rotating magnetized gaseous bodies. For example, convection may occur due to 
an outwardly decreasing entropy in radiatively inefficient magnetized accretion flow (Narayan 
et al. 2002) or in radiation-dominated accretion disks (see, e.g., Agol et al. (2001). If the 
angular velocity depends on the vertical coordinate, then stratification can be the reason of 
the baroclinic convection associated with non-parallel temperature and pressure gradients 
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(Urpin & Brandenburg 1998, Urpin 2003). Stratification can also lead to a wide variety 
of strong non-axisymmetric instabilities in accretion disks (Keppens, Casse, & Goedbloed 
2002). 

Many previous stability analyses have adopted the Boussinesq approximation, and have 
therefore neglected the effect of compressibility. This is valid provided that the magnetic 
field strength is essentially subthermal, and the sound speed is much greater than the Alfven 
velocity, c s ^> ca- However, this can often be not the case in both astrophysical objects and 
many numerical simulations. The first attempt to consider the effect of compressibility on 
the magnetorotational instability was undertaken by Blaes & Balbus (1994) in the context 
of astrophysical disks. The authors consider a very simplified case of the wavevector parallel 
to the rotation axis and vanishing radial magnetic field. As a result, the most interesting 
physics has been lost in their study since only the standard magnetorotational instability 
operates in this simple geometry. For instance, the authors concluded that the azimuthal 
field does not affect the instability criterion in all orders in B^. This conclusion applies, 
however, only to the particular case considered by Blaes & Balbus (1994) and is incorrect in 
a more general magnetic configuration. A detailed study of the influence of B^ on rotation 
induced instabilities is given by Pessah & Psaltis (2005). These authors show that terms in 
the Lorentz force proportional to B^/s where s is the cylindrical radius can be dynamically 
important when superthermal magnetic fields are considered. The curvature of toroidal field 
lines plays a fundamental role in the local stability of strongly magnetized compressible 
plasmas. In particular, the combined influence of compressibility and the curvature term 
can give rise to two new instabilities. 

A more refined analyses of instability of a differentially rotating magnetized gas was 
undertaken by Blokland et al. (2005) and van der Swaluw, Blokland, & Keppens (2005) 
who considered radially stratified disks in the cylindric limit. In particular, Blokland et 
al. (2005) investigate the influence of a toroidal field on the growth rate and frequency of 
eigenmodes and find that it leads to overstability (complex eigenvalue) and reduction of the 
growth rate. The overstable character of the magnetorotational instability increases as the 
angular velocity decreases. Van der Swaluw, Blokland, & Keppens (2005) study the interplay 
between different instabilities in accretion disks and argue that the growth rate of convection 
can be essentially increased due to magnetorotational effects. Note that both these studies 
treat stability in the magnetic field with a vanishing radial component. 

In this paper, we show that a number of instabilities may occur in a compressible dif- 
ferentially rotating magnetized gas which are different from the standard magnetorotational 
instability. These instabilities appear if the magnetic field has non-vanishing radial and 
azimuthal components and can manifest themselves even if stratification is neglected. The 
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considered instabilities can be oscillatory or non-oscillatory and arise even if the magnetic 
field is strong enough to suppress the magnetorotational instability. Our study generalizes 
the paper by Bonanno & Urpin (2005) where shear-driven magnetic instabilities were con- 
sidered in a compressible flow for the particular case when the wavevector is perpendicular 
to the magnetic field and the standard magnetorotational instability does not operate. 

The paper is organized as follows. In section 2, we discuss the main equations and derive 
the dispersion relation for a compressible differentially rotating magnetized gas. Stability 
properties of modes are considered in section 3, and our numerical calculations of the growth 
rate are presented in section 4. Finally, we discuss our conclusions in section 5. 



2. Basic equations and dispersion relation 

We work in cylindrical coordinates (s, ip, z) with the unit vectors (e s , e v , e z ). The 
equations of compressible MHD read 

V» 1 - 

v + (v-V)v = - +g + —(V x B) x B, (1) 

p Anp 

p + V-(pv) = 0, (2) 

p + v ■ Vp + 7pV • v = 0, (3) 

B-Vx(ijxB) = 0, (4) 

V-B = 0. (5) 

Our notation is as follows: p and v are the gas density and velocity, respectively; B is the 
magnetic field; p is the pressure and 7 is the adiabatic index. 

The basic state on which the stability analysis is performed is assumed to be quasi- 
stationary with the angular velocity Q = Q(s) and B 7^ 0. We assume that gas is in 
hydrostatic equilibrium in the basic state, then 

— = D + -!-(V xB)xB , D = g + fl 2 s. (6) 
p A-np 

We consider magnetic configurations where both the radial and azimuthal field components 
are presented. The presence of a radial magnetic field and differential rotation in the basic 
state can lead to the development of the azimuthal field. Nevertheless, the basic state can 
be considered sometimes as quasi-stationary despite the development of the toroidal field. 
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For example, if the magnetic Reynolds number is large, one can obtain from Eq. (4) that 
the azimuthal field grows approximately linearly with time, 

B v (t) = B v (0) + stfB s t, (7) 

where Q' = dVt/ds, and -6^(0) is the azimuthal field at t — 0. As long as the second term on 
the r.h.s. is small compared to the first one or, in other words, 

stretching of the azimuthal field does not affect significantly the basic state; r v is the char- 
acteristic timescale of generation of B^. As a result, the basic state can be treated as 
quasi-stationary during the time t < r^. If B v (0)/B s ^> 1, then steady-state can be main- 
tained during a relatively long time before the generated azimuthal field begins to influence 
the basic state. 

We consider stability of axisymmetric short wavelength perturbations. Small perturba- 
tions will be indicated by subscript 1, while unperturbed quantities will have no subscript. 
The linearized momentum equation is 

vi + (vi ■ V)v + (v • V)vi = 



P 

Pi r ""~ 



P 

where 



^ - —(V xB)xB 
p Anp 



+ L, (9) 



[(VxBi)xB+(VxB)xBi], 



Taking into account that the unperturbed motion is rotation {v = sfle v ) and using Eq. (6), 
we can transform this equation into 

b 1 + 2nxv 1 +e v stfv ls = -^ + ^D + L. (10) 

P P 

The coefficients of Eq. (10) are not constant, but they are slowly varying functions of coor- 
dinates compared to perturbations. Therefore, one can expand coefficients in Taylor series 
around some radial point and retain only the zeroth order terms in this expansion. This is 
a good approximation as long as the wavelength of perturbations is small compared to the 
lengthscales of unperturbed quantities. Then, one can perform a local stability analysis of 
Eq. (10) around this radius but with constant coefficient. Since coefficients are constant, 
one can consider perturbations with the spacetime dependence oc exp(at — ik ■ r). where 
k = (k s , 0, k z ) is the wavevector. 
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The term proportional to p±/p on the r.h.s. of Eq. (10) is typically small in the short 
wavelength approximation. Indeed, we can estimate pi ~ c 2 s pi where c s is the sound speed. 
Then, the pressure term in Eq. (10) is kiC 2 (pi/ p) where ki is the cylindrical component of 
the wavevector. Therefore, the ratio of the first two terms on the r.h.s. of Eq. (10) can 
be estimated as kiC 2 s /Di where D = g + Q 2 s. In many cases of interest, the centrifugal 
force is weak compared to gravity, and we have Di ~ gi. This estimate applies in stellar 
radiative zones, solar tachocline, protostellar disks, galaxies, etc. However, this estimate 
can be incorrect in accretion disks with substantially superthermal magnetic fields (Pessah 
& Psaltis 2005), but we do not address these objects in the present paper. Since <?J g% is 
the pressure scaleheight H^, we obtained that the ratio of the first and second terms on the 
r.h.s. of Eq. (10) is ~ k{Hi. The inequality kiHi 3> 1 is exactly equivalent to the condition of 
applicability of the local approximation.. The terms of the order of (/cjifj) -1 cannot be taken 
into account consistently in the local approximation because inhomogeneities of the basic 
state produce corrections of the same order of magnitude. Therefore, we have to neglect the 
term proportional to D in what follows. 

Consider in more detail the Lorentz force in Eq. (10). For the sake of simplicity, we 
assume that the radial and azimuthal components the magnetic field in the basic state do 
not depend on the vertical coordinate, dB^/dz = dB s /dz = 0, and that dB z /ds = 0. Then, 
the electric current in the basic state has only the vertical component, and the Lorentz force 
in Eq. (10) can be represented as 



4npL = e s B z 

s 2 ds \ v 
oz s 



dB u dB lz 
dz ds 
B ls d(sB^) 



B u 



dB 



ds 



ds 



-e. 



dB u 



' ds 
dB lz 



dz 



+B a 



ds 



dB 



lip 



dz 



The last terms in the s- and (/2-components of this equation are of the order of (ks)^ 1 
compared to the second terms and should be neglected in a short wavelength approximation 
that does not allow to treat consistently corrections of such order because a slow dependence 
of the basic state on coordinates (which is neglected in a local analysis) produces corrections 
of the same order. Then, the Lorentz force with the accuracy in terms ~ X/s, where A is the 
wavelength of perturbations, is given by 



1 



4np 



iBx(kxB 1 ) + ^—^ 1 



ds 



(12) 
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The second term in this expression describes the effect caused by the presence of electric 
currents in the basic state. This term can be important even in a local approximation if the 
toroidal field is sufficiently strong and is not proportional to s _1 . 

Following similar transformations with the remaining equations in the system given 
by Eq. (l)-(5), we arrived to the linearized equations needed for stability analysis. These 
equations read with accuracy in the lowest order in X/s 



av\ + 2Q x v\ + e v sQ!vi s = 

P 

' (fcxB 1 )xE + ^^ ) (13) 



4np 4irp s ds ' 

ap 1 -ip(k-v 1 ) = 0, (14) 

api — i~fp(k ■ Vi) = 0, (15) 

(7 B\ = e v sVt'B ls - i(B ■ k)^ + iB{k ■ v^, (16) 

k-B 1 = 0. (17) 

Assuming B v < (ks)B s , Eqs. (13)-(17) may be combined after some algebra into a sixth-order 
dispersion relation, 

o" 6 + a4<7 4 + a3<r 3 + a 2 <T 2 + a\<j + a = 0. (18) 

The coefficients of this equation can be expressed in terms of characteristic frequences. We 
have with accuracy in the lowest order in ~ [ks) v 

a 4 — (Jq + u 2 A + f^, a 3 = + 2iilil Aip u; A , 
a 2 = i0 2 A (u 2 s + Uq) + fi^/iWo + uj A k s c As ) 

—AVt 2 UJAk z CAz + i^lA^A c A<fik s , 

ai = uj% n (u A + p^l) + 2iVLVL Aip oj A [p^l + 

UA(k s c As - k z c Az )}, 
a = (uj 2 u 2 a + itt A(p uJ A CAsCA V k 2 )(uJ 2 A + pQ 2 sh ), 

where \i = k 2 /k 2 . The characteristic frequencies are given by 

n 2 e = 4n 2 + n 2 sh , n 2 sh = 2 S nn' , 

2 2,2 2 2r2 2 2;2 

uj =uj s + u m , uj s = c s k , uj m = c A k , 
lu a = k ■ c A , c A = B/ , 

^bq = k 2 c Alfi CA S sn', n Av = ~t^( 1+ ^ > 
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where c s = \fjpfp is the sound speed. Eq. (18) describes six modes that can generally exist 
in a compressible rotating magnetized gas. 

Before to solve Eq. (18), we consider several particular cases and compare our dispersion 
relation with those derived previously by other authors. In the incompressible limit when the 
sound speed c s and, hence, the sound frequency uj s are very large, we obtain from Eq. (18) 
the dispersion relation for the standard magnetorotational instability, 

a 4 + a 2 (2u 2 A + + cu 2 A (cu 2 A + ^ 2 sh ) = (19) 

(Velikhov 1959). This equation allows unstable solutions if 

fin 2 ah < -u 2 A , (20) 

that represents the criterion of the magnetorotational instability in the incompressible limit. 
It is seen that the magnetorotational instability can be entirely suppressed if the magnetic 
field is sufficiently strong (see, e.g., Urpin 1996, Kitchatinov & Rudiger 1997, Balbus & 
Hawley 1998). 

In the particular case k s = B s = and d(sB v )/ds = (no currents in the basic state), 
we obtain from Eq. (18) the dispersion relation derived by Blaes & Balbus (1994) in their 
analysis of the effect of compressibility, 

a 6 + q^ + q 2 (T 2 + q = 0, (21) 

where 

q A = k 2 (2c 2 Az + c 2 M + c 2 s ) + n 2 e , 
q 2 = k 2 [k 2 c 2 Az (c 2 A + 2c 2 s ) + Q 2 e (c 2 s + c 2 M ) + c 2 Az Q 2 sh ], 

qo = k*c 2 s c 2 Az (k 2 c 2 Az + n 2 h ). 

From Eq. (21), Blaes & Balbus (1994) deduced that instability of compressible and incom- 
pressible gas is determined by the same criterion (20). 

Note that Eq. (18) in the case B s = differs from the dispersion relation derived by 
Pessah & Psaltis (2005) (see Eq. (25) of their paper). These authors consider the influence of 
the "curvature terms" on stability of differential rotation in accretion disks. Likely, the dif- 
ference is caused mainly by different assumptions regarding the basic state. Pessah & Psaltis 
(2005) derive the dispersion relation for the midplane of accretion disks and a substantially 
superthermal magnetic field. As it was mentioned above, our study is not addressed such 
objects. 
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Eq. (18) is very complicated for analysis, and we consider in this paper only the case 
when the effect of electric currents can be neglected. Certainly, this is valid if B v oc 1/s, but 
the contribution of currents can be small in some other well. Assuming sQ' ~ Q, 

we can estimate that the terms caused by the electric current (oc fi^) in the coefficients a 3 
and a2 provide only small corrections of the order of A/s to uj B q and Wqu\, respectively, and 
can be neglected in these coefficients. The effect of currents in a\ is unimportant is 

B% « (ks)B 2 . (22) 

Finally, the influence of currents in the coefficient a is negligible if 

4 < ks, (23) 

that is fulfilled usually in astrophysical conditions. Conditions (22) and (23) determine the 
domain of parameters where the effect of electric currents on stability properties can be 
neglected. For this reason, we consider in the present paper stability in the domain given by 
Eqs. (22) and (23) and will suppose Qa<p = in Eq. (18). 

It was shown by Bonanno & Urpin (2006) that Eq. (18) at Qa<p — has unstable solutions 
even if the criteria of the magnetorotational and Rayleigh instability are not satisfied. The 
authors consider stability of perturbations with perpendicular k and B. The dispersion 
relation reads in this case 

a 5 + a\uj 2 Q + Q 2 e ) + a 2 cu 3 Bn + a^W + ^"in = 0. (24) 

— * — * 

As it is seen from Eq. (19), perturbations with k ■ B = are not subject to the magnetoro- 
tational instability in the incompressible limit. Nevertheless, as it was shown by Bonanno & 
Urpin (2006), Eq. (24) can have unstable solutions if uj B q ^ 0. Depending on parameters, 
the roots of Eq. (24) can vary within a wide range, but the growth rate of unstable modes 
is always proportional to a power of ubq (if the magnetorotational and Rayleigh instability 
do not occur in a flow). 

The presence of terms oc ujbq in Eq. (18) is always crucial for instability. To illustrate 
this point, we consider the solution of dispersion relation (18) in the case of a non-vanishing 
Alfven frequency, uja ^ 0. For the sake of simplicity, we assume that Vt 2 = (Q oc s~ 2 ) and 
consider the case of a cold gas with small c s . Since the magnetorotational instability does 
not occur if ca 3> c s , the flow is not subject to any known rot at ion- induced instabilities. 
Dispersion relation (18) reads in this case 



a 5 + a 3 (u 2 m + u 2 A ) + <7 2 c4 n + a{u 2 A u 2 m 
-4Q 2 u A k z c A z} + Ua^bu = 0- 



(25) 
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We consider stability of a special type of perturbation with the wavevector, satisfying the 
condition k 2 = AhVL 2 /c A where h = k z CA z /^A- If sQ/ca 1, our choice of k is compatible 
with a short wavelength approximation since it implies ks >> 1. For such perturbations, we 
have from Eq. (25) 

<r V + u4 + c4) + (a 2 + cv 2 A )u 3 BQ = 0. (26) 

If the radial component of the magnetic field is small, B s <C B z ~ 5^, then the frequency 
uosn can be considered as a small perturbation in Eq. (26), and we can solve this equation by 
making use of the perturbation procedure. Two roots are non- vanishing even in the zeroth 
order, and they can be obtained approximately by neglecting ubq in Eq. (26). Then, we 
have 



<7l,2~±^^ + ^. (27) 

The growth rate of oscillatory modes is small compared to the frequency in this case, and 
we neglect it in the expression for a± t 2- Three other roots are proportional u Bn , and their 
frequencies are small compared to u m and u A . Therefore, calculating these roots, we can 
neglect a 2 compared to uo 2 m and uj 2 a in both brackets in Eq. (26). Then, we obtain 

3 ' 4 ' 5 ~ LU 2 + LU 2 ' ^ ' 

U A + u m 

Since oja ~ uj m , the growth rate of unstable modes can be estimated as Re a ~ |w_bq|. Taking 
into account that k 2 ~ f2 2 /c 2 A , we obtain 

/ B \ 1/3 

Recr-fl y-^-j (29) 

(we assume that sQ! ~ Q). Note that similar solutions with Re a ~ \oubci\ exist also 
in flows with non- zero sound speed if c s satisfies the condition c s <^ min(c J 4 s , Vt/k) where 
c As = Bsj^fWp. 



3. Criteria of instability 

The conditions under which Eq. (18) has unstable solutions with a positive real part can 
be obtained by making use of the Rooth-Hurwitz theorem (see Henrici 1977, Aleksandrov, 
Kolmogorov, & Laurentiev 1985, Toth, Szili, & Zachar 1998). The Rooth-Hurwitz criteria 
for a sixth-order equation are rather complex in the general case, but they are very much 
simplified for Eq. (18) where the coefficient before a 5 is vanishing. In accordance with these 
criteria, Eq. (18) has unstable solutions if one of the following inequalities is fulfilled 



a 3 > , 



(30) 
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a\ > 

01(0304 — ai) — 02^3 < , 
a\(a z a A — a x ) + 03(0003 — Oia 2 ) < , 

a < . 

One of conditions (30) and (31) is always fulfilled if a 3 7^ 0, or 



(31) 
(32) 
(33) 
(34) 




(35) 



Apart from differential rotation, this criterion requires non- vanishing radial and toroidal com- 
ponents of the magnetic field. Condition (35) can generally be satisfied for both the inward 
and outward decreasing angular velocity in contrast to the criterion of the magnetorotational 
instability (19) that is fulfilled only if ft' < 0. 

Since Eq. (18) has always unstable solutions if ujbq 7^ 0, remaining criteria (32)-(35) 
can be treated only in the case ubq = 0. In this case, however, we have a\ = 03 = 0, 
and neither Eq. (32) nor Eq. (33) is fulfilled. Therefore, these criteria do not provide any 
new conditions of instability. Eq. (34) yields the standard criterion of the magnetorotational 
instability (19) which does not overlap condition (35). For instance, the magnetorotational 
instability can occur even if B s = or = 0, when the magnetic shear-driven instability 
(35) does not operate. On the contrary, the shear-driven instability can arise under conditions 
when the magnetorotational instability does not occur. For example, as it was mentioned, 
the magnetorotational instability is suppressed by a sufficiently strong magnetic field (see 
criterion (19)) but instability (35) does not exhibit such suppression and may occur even if 
the magnetic field is relatively strong (see Eq. (23)). 

It should be emphasized that differential rotation always generates the azimuthal mag- 
netic field from the radial one and, in fact, the necessary and sufficient condition of instability 
(35) is a non-vanishing radial magnetic field. Therefore, we have to conclude that all differ- 
entially rotating gaseous flows with a non-vanishing radial magnetic field are unstable, and 
this conclusion does not depend on the sign of B s or ft'. 

Note that there is no one-to-one correspondence between the roots and criteria. In some 
cases, few modes are unstable if only one criterion is fulfilled and, on the contrary, only one 
mode can be unstable if few criteria are satisfied. Generally, Eq. (18) describes two Alfven 
waves, two slow and two fast magnetosonic waves, however, it can be difficult to discriminate 
between these modes if ft ~ oj s ~ uja- To demonstrate which modes are unstable if condition 
(35) is satisfied we consider a very slow rotating gas with u s ~ uja ft ~ sft'. In this case, 
we can restrict the solution of Eq. (18) only by terms of the zeros and first orders in ft and 
ft'. Then, Eq. (18) yields for such slow rotation 



(a 2 + iu 2 A )[(a 4 + a 2 u 2 Q + lu 2 lu 2 a ) + aco 3 Bn \ = 0. 



(36) 
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The term oc u% n describes a small linear corrections associated to differential rotation. The 
magnetorotational instability does not appear in the linear approximation in Q since it is 
caused by quadratic in the angular velocity terms. A couple of roots of Eq. (36) corresponds 
to to the Alfven waves, a = ±iouA, which are stable in the linear approximation even if 
criterion (35) is fulfilled. Note, however, that Alfven waves can be unstable if non-linear in 
Q terms are taken into account. Fast and slow magnetosonic modes satisfy the dispersion 
relation of the fourth order, 

4,2 2,2 2 3 /o7\ 

a + a u + u s u A = -cru Bn . (37) 

This equation can be solved by making use of the perturbation procedure since is small. 

We can represent a as cr + Act where do is the solution of the dispersion equation at Q = 
and Act is small correction caused by slow rotation. Then, it can be obtained from Eq. (37) 
that 

-o 2 = ± ^--M, (38) 

and 

2Aa = =p4n("o - ^M)" 1/2 - (39) 

The imaginary part of a is determined by o"o whereas the real part which is responsible for 
instability is determined by Act. The upper and low sign in Eq. (39) corresponds to the 
slow and fast magnetosonic modes, respectively. Only these modes can be unstable under 
condition (35) in the linear in Q approximation. If < 0, then the slow mode is unstable 
and, on the contrary, if u>% n > then the fast mode is unstable. Note, however, that this 
simple picture can be changed if rotation is fast, fl ~ uj s ~ lo\. The Alfven mode can become 
unstable in this case as well. 

Our conclusion regarding instability caused by a combined effect of compressibility and 
differential rotation in a magnetized gas appears to be in contrast to the result obtained 
in the review by Balbus & Hawley (1998). These authors claim that in the presence of 
Keplerian differential rotation, the fast magnetosonic mode remains stable, and only slow 
magnetosonic mode can be unstable. In fact, however, there is no contradiction between these 
conclusions because Balbus & Hawley (1998) consider the basic magnetic configurations with 
only vertical and azimuthal magnetic field components, but the radial field is vanishing in 
their model. As it follows from our analysis, ubq = in this case, and the instability 
associated to compressibility should not occur in the model considered by Balbus & Hawley 
(1998). 
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4. Numerical results 

To calculate the growth rate of instability, it is convenient to introduce dimensionless 
quantities 

a_ = 4^ ur[l_ = a4j_ 

fi e ' fi2 ' h fi2 x 2 ' M ^2 x 2 > 

_ (fc • £) 2 , _ C^CAs , _ A^CAz 
k 2 B 2 c A uja 

(we assume Vt 2 > 0). Note that the sign of tt' is determined by the parameter q: £1' is 
positive if q < 1 and negative if q > 1. Then, Eq. (18) transforms into 

r 6 + r 4 [i + / s x 2 + f mX \\ + e)] + r 3 ^-5f m x 2 

+r¥[e/ m (2/ s x 2 + f m x 2 + 1 - - hq) 
1 — q 

+Kfs + fm)} + r — — (/x + ef m x 2 )5f m x 2 

+ef s f m x 4 [ef m x 2 + /i(l - g)] = 0. (40) 

The dependence on the wavelength is characterized by the parameter x 2 in this equation. 
Eq. (40) was solved numerically for different values of the parameters by computing the 
eigenvalues of the matrix whose characteristic polynomial is given by Eq. (40) (see Press et 
al. 1992 for details). 

In Fig. 1, we plot the dependence of the real part of T on f m for the flow with q = 4 
(Q oc s~ 3 / 2 ) and /i — e — 0.3, 5 = 0.2, h = 0.5 and x 2 = 10. Solid and dashed line correspond 
to complex and real roots, respectively. We choose the value f s = 10 for Fig. 1. Therefore, 
fs f m for the considered range of f m and, hence, c s ^> C4 that corresponds to the case 
when compressibility is almost negligible. The profile of Q with q = 4 is unstable to the 
magnetorotational instability, and our calculations show that this instability can actually 
occur if f m is not large. Criterion of the magnetorotational instability (19) can be rewritten 
in terms of the dimensionless parameters as 

f m < 3fi/ex 2 . (41) 

For given values of fx, e, and x 2 , we obtain that this inequality satisfies in the region where 
f m < 0.3. It is seen from Fig. 1 that there is a real root with large Re T (~ 0.1 — 0.4) that 
can be identified with the magnetorotational instability (which is always non-oscillatory in 
the incompressible limit). Apart from this root, however, there exist another unstable root 
which is much smaller for the chosen values of parameters. The new root is positive in the 
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region f m > 0.3 where the magnetorotational instability is completely suppressed. The new 
instability exists even in the region of very large f rn in complete agreement with criterion 
(32). The important point is that this unstable root is caused by compressibility and does 
not appear in the Boussinesq approximation. The growth rate of the new instability is 
relatively small (Re Y ~ 0.01) for the chosen parameters because c s ^> ca, and departures 
from the incompressible limit are small. The situation, however, can change qualitatively if 
the magnetic field is stronger, and we are beyond the necessary condition of the Boussinesq 
approximation, ca <C c s . 

In Fig. 2, we show the same dependence as in Fig. 1 but for the case f s = 0.3 (other 
parameters are same) when compressibility plays much more important role. The Boussinesq 
approximation does not apply in this case except the region of small f m <C f s . Nevertheless, 
the magnetorotational instability can still occur in the region f m < 0.3, and its growth rate 
is in a good agreement with the incompressible case. In the region of a strong field, however, 
there exists a new instability caused by shear and the magnetic field. Compressibility is 
crucial for this new instability, and it cannot be obtained in the incompressible limit. The 
growth rate of the newly found instability is rather high, Re T ~ 0.1 — 0.2, in the region 
f m > 0.3 where the magnetorotational instability is suppressed. At a fixed wavelength of 
perturbations, the growth rate increases slightly if the magnetic field becomes stronger. The 
new instability is not suppressed even by a very strong magnetic field. Calculations show 
that further decrease of f s (that is equivalent to increase of compressibility) leads to a higher 
growth rate that can be even as high as Re T ~ 0.5 in some cases. 

The number of unstable roots and the growth rate are sensitive to a relative orientation of 
the azimuthal and radial magnetic fields. This orientation is characterized by the parameter 
5, and we plot the growth rate versus f m for 5 = —0.2 in Fig. 3. In contrast to Fig. 2, there 
are two couples of unstable complex conjugate roots for 5 = —0.2. As it was mentioned in 
Sec. 3, apart from magnetosonic waves the Alfven mode can be unstable as well if rotation is 
sufficiently fast. One couple has a larger growth rate (Re T ~ 0.1) whereas another couple 
grows substantially slower (Re T ~ 0.05). As usual, compressible instabilities are most 
efficient in the region where the magnetorotational instability does not occur (f m > 0.3). 
For 5 = —0.2, however, the compressible instabilities can arise even at f m < 0.3, where they 
work in parallel to the magnetorotational instability. 

The above three figures illustrate the behavior of roots in the case when Q' < and, 
hence, the magnetorotational instability can operate if the magnetic field satisfies Eq. (41). 
The newly found instability, however, can arise in the case Q' > as well. In Fig. 4, we plot 
the real parts of roots for q = 0.5 that corresponds to f2 oc s 2 (other parameters are same 
as in Fig. 1). It is seen from the figure that there is no large root that could correspond 
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Fig. 1. — The dependence of the real part of r on f m for /x = e = 0.3, 5 = 0.3, h = 0.5, 
x 2 = 10, and / s = 10. The angular velosity profile is given by il oc s~ 3//2 . 




Fig. 2. — Same as in Fig. 1 but for f s = 0.3. Solid and dashed lines correspond to complex 
and real roots, respectively. 
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to the magnetorotational instability in the region of low magnetic fields. Instead, there are 
three couples of complex conjugate roots, and one of these couples has positive Re T. These 
roots represent a relatively weak shear-driven instability that operates at any field strength. 
Obviously, the instability is weak (Re T ~ 0.01) because c s ^> c\ (or f s ^> f m ) for the chosen 
parameters. However, the growth rate becomes larger as the magnetic field increases and 
departures grow from the incompressible limit. 

In Fig. 5, we show the growth rate for the case when the magnetorotational instability 
still does not occur (q = 0.5 and Q oc s 2 ) but the effect of compressibility is more significant 
than in Fig. 4 since f s = 0.3. Again, all modes are complex, and two couples of complex 
roots are unstable. Due to compressibility, the growth rate of the most unstable modes can 
reach rather high values ~ 0.04 — 0.06f2 e . Since the epicyclic frequency is large, VL 2 = 8Q 2 , we 
estimate the growth rate as 0.1 — 0.2Q, that implies that the growth time of instability is of 
the order of the rotation period, P = 2n/Q, in this case. Therefore, the considered instability 
turns out to be very rapid even in those cases when the magnetorotational instability cannot 
operate. Note that the newly found instability may occur even if the magnetic field is very 
strong and the magnetic pressure is greater than the gas pressure (or f m > f s ). 

Generally, a dependence of the growth rate on the parameters \i and e is relatively 
weak whereas the value h plays more important role for instability. In Fig. 6, this point is 
illustrated by calculations of Re T for h = 2 that is larger than in the previous figures. In 
fact, the parameter h characterizes the direction of k with respect to the magnetic field. The 
value h = 2 implies that k z /k s = —2B S /B Z . Calculations have been performed for f s = 0.3 
to compare with the previous figure. A higher value of h does not change complexity of the 
modes, and all roots are complex. Generally, an increase in h can lead to a larger number of 
unstable modes and a higher growth rate as it is seen from the figure. The dependence of Re 
T on the field strength becomes non-monotonic with maximum at f m ~ f s . The maximum 
growth rate in Fig. 6 is approximately 3 times greater than in the case h = 0.5. 

In Fig. 7, we plot the dependence of the real and imaginary parts of T on x 2 for the 
case of radially decreasing Vt(s) (q = 4) and for v = e = 0.3, 5 = 0.2, h = o.5, f s = 0.3, and 
f m = 0.2. The unperturbed field strength is not sufficient to suppress the magnetorotational 
instability. Solid lines show the growth rate and frequency for the complex roots, and dashed 
lines for the real roots. The real roots correspond to the magnetorotational instability which, 
in accordance with condition (41), can occur only if x 2 < 15 for the given parameters. In 
the region of shorter wavelengths (x 2 > 15), however, the magnetorotational instability is 
forbidden but the newly found instability can be efficient. In this region, there are two 
pairs of unstable complex conjugate roots and a pair of stable roots. The new instability is 
oscillatory, and arises despite the magnetic field is rather strong (the thermal and magnetic 
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Fig. 4. — Same as in Fig. 1 but for the case q = 0.5 when the magnetorotational instability 
does not occur. 
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Same as in Fig. 4 but for / s = 0.3. 
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energy are comparable). Even in such a strong field, the growth time of instability is rather 
short, ~ 3P. Note that the growth rate is weakly dependent on the wavelength for x 2 > 15. 

In Fig. 8, we plot the same dependence as in Fig. 7 but for the flow with q = 0.5 which 
is stable to the magnetorotational instability at any magnetic field. Only complex modes 
exist in such a flow and, as a result, instability is oscillatory. In this case, Re T is a bit 
smaller than for a flow with the Keplerian angular velocity (q = 4) but, because of a high 
epicyclic frequency Q e = 2y/2D,, the growth time is rather short, ~ P. As in the previous 
figure, the growth rate is approximately same for all short wavelength perturbations. 

Finally, we plot in Fig. 9 the dependence of the growth rate on x 2 for rapidly rotating 
gas with f m = f s = 7 x 10 -3 . This figure illustrates the type of solutions that can exist 
in a relatively cold plasma. For such solutions, the growth rate is rather high, ~ wbci- 
Remarkably, that this high growth rate is achieved for perturbations with k 2 ~ Q 2 /c 2 A or, 
using the dimensionless parameters, x 2 ~ 1/ f m . The calculated growth rate is in a reasonable 
agreement with estimate (26) since Bg/B^ 0.1 in the considered case. The growth rate 
is smaller than it is given by Eq. (29) for very large x 2 . Note that the roots (27)-(28) were 
calculated under assumption ca 3> c s whereas c s = ca in the numerical solution. Therefore, 
even the number of modes is different in Eqs. (24)-(25) and in Fig. 9. However, estimate of 
the growth rate (29) is still valid despite this difference. 



5. Discussion 

We have considered stability of differentially rotating magnetized gas taking account of 
the effect of compressibility and assuming that the magnetic field has non-vanishing radial 
and azimuthal components. In our stability analysis, we assume that the basic state is quasi- 
stationary. This assumption can be fulfilled in many cases of astrophysical interest despite 
the development of the azimuthal field from the radial one due to differential rotation. For 
instance, if the magnetic Reynolds number is large then the timescale of generation of the 
toroidal field is ~ (see Eq. (8)) in the case B v (0) > B s . Obviously, the basic state can be 
considered as quasi-stationary if the growth time of instability is shorter than t v . As it is 
seen from Eq. (29), the growth rate of instability can be roughly estimated as Q(B S /B (p ) 1 / 3 
if sQ' ~ Q. Then, the condition of quasi- st at ionarity reads 

B V ^>B S . (42) 

The physics of the instability is very simple. To clarify its nature, we consider the 
simplest case of a cold gas with ca c s . Assume that the background magnetic field has 
only radial B s and azimuthal B v component, and the perturbation is a plane wave in the 
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Fig. 7. — The dependence of the growth rate and frequency on x 2 = (ks) 2 for radially 
decreasing Q(s) (q = 4) and for ji — e — 0.3, 5 = 0.2, h = 0.5, f s = 0.3 and f m = 0.2. 
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Fig. 8. — Same as in Fig. 7 but for q = 0.5. 
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z-direction, such as k = (0,0, k z ). The wave produces the density perturbation p\ that 
is also an oscillatory function of z. Since the magnetic field is assumed to be frozen into 
the gas, compression/decompression should perturb the magnetic field. The perturbation 
of the radial component is simply related to the density perturbation, B ls = B s (pi/p). 
The azimuthal perturbation is determined by two processes and can be written as B\ v = 
B^+Bty, where B$ = B ip (pi/p) is the perturbation caused by compression/decompression 
in the wave and B[ s ^ is caused by stretching of the azimuthal field from the perturbation 
of the radial component due to differential rotation. The rate of stretching is sQ'B ls and, 
hence, b[ s ^ = sQ' J B ls dt'. Then, the perturbation of the azimuthal field is given by 



B\ v = B v — + sQ 
P 



'B s J ^dt'. (43) 



Note that this equation is exactly equivalent to the (^-component of Eq. (16) if we assume 
k ■ B = and substitute i(k ■ V\) = p\j p as it follows from the continuity equation (14). 
The perturbation of the field leads to the Lorentz force which, in its turn, generates the 
vertical velocity. The Lorentz force is equal to (V x B ± ) x B/An = ik(B ■ B^/An. Then, the 
perturbation of the vertical velocity is 

In its turn, vertical motions should generate perturbations of the density since the gas is 
compressible. We have from the continuity condition 

^ = ik z J v lz dt'. (45) 

Combining Eqs. (43)-(45), we obtain the expression governing density perturbations 

7 = / *7 *" + sa ' B ' B * I p *"') • (46) 

If we assume that p\ oc exp(crt), then the growth rate is determined by the condition 

a 3 + auj 2 + = 0. (47) 

— * 

This equation describes the same magnetic modes as Eq. (24) does in the case k = k z e z . 

— * 

Indeed, in the case k = (0, 0, k z ) (or p — 1), Eq. (24) can be transformed into 

(a 2 + Q 2 e )(a 3 + au; 2 +u; 3 Bn ) = 0, (48) 



such as the modes affected by the magnetic field are described by Eq. (47). It is seen from 
this simple explanation, that shear and compressibility play the major role in the instability. 
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In the incompressible limit that applies if c s ^> ca, differentially rotating flows can 
be subject to the well-known magnetorotational instability (Velikhov f 959, Chandrasekhar 
1960) that arises if the angular velocity decreases with the cylindrical radius. It turns out, 
however, that compressibility alters drastically the stability properties of magnetized rotating 
flows. Apart from instabilities caused by "the curvature terms" (see Pessah & Psaltis 2005), 
a number of new instabilities can occur in the presence of a radial magnetic field. The 
properties of these instabilities are very much different from those of the magnetorotational 
instability. The necessary condition of these new instabilities is given by Eq. (35) and reads 
BsB^Q! 7^ 0. Since differential rotation always generates the azimuthal magnetic field from 
the radial one, the necessary and sufficient condition of instability is the presence of a radial 
field in differentially rotating flows. Obviously, this condition can widely be satisfied in 
astrophysical bodies. Note that, contrary to the magnetorotational instability that occurs 
only if the angular velocity decreases with the cylindrical radius, the newly found instabilities 
may occur at any sign of Q! . 

Likely, the most important difference between the magnetorotational and newly found 
instabilities is associated with the dependence on the magnetic field strength. A sufficiently 
strong magnetic field, satisfying the inequality (k ■ B) 2 > 87rpsQ\Q'\(k 2 /k 2 )\, suppresses 
completely the standard magnetorotational instability. On the contrary, the instability found 
in our paper cannot be suppressed even in relatively strong magnetic fields (but satisfying 
inequality (23))as it is seen from the criterion (35). Often, the growth rate reaches some 
saturation in very strong magnetic field, and this saturated value can be rather high (Re 
r~ 0.1 -0.5). 

The growth rate of the newly found instabilities is rather high and can even be com- 
parable to ~ Q e for some values of parameters. Generally, the growth rate depends on 
compressibility, being smaller for a low compressibility. The incompressible limit (Boussi- 
nesq approximation) corresponds to c s ^> ca, and the considered instabilities are inefficient 
in this limit because of a low growth rate. However, in the case of a strong field with 
Ca ~ c s when the Boussinesq approximation does not apply, the instability can be much 
more efficient than the magnetorotational instability. 
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